Roughness at the depinning threshold for a long-range elastic string 
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In this paper, we compute the roughness exponent (^ of a long-range elastic string, at the depinning 
threshold, in a random medium with high precision, using a numerical method which exploits 
the analytic structure of the problem ('no-passing' theorem), but avoids direct simulation of the 
evolution equations. This roughness exponent has recently been studied by simulations, functional 
renormalization group calculations, and by experiments (fracture of solids, liquid meniscus in *He). 
Our result C, = 0.390 ± 0.002 is significantly larger than what was stated in previous simulations, 
which were consistent with a one-loop renormalization group calculation. The data are furthermore 
incompatible with the experimental results for crack propagation in solids and for a *He contact 
line on a rough substrate. This implies that the experiments cannot be described by pure harmonic 
long-range elasticity in the quasi-static limit. 
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The statics and dynamics of elastic manifolds in ran- 
dom media govern the physics of a variety of systems, 
ranging from vortices in type-II superconductors H and 
charge density waves |2| to interfaces in disordered mag- 
nets 111, contact lines of liquid menisci on a rough sub- 
strate^H] and to the propagation of cracks in solids p] . 

In most cases, the restoring elastic forces acting on 
a point on the manifold are local i.e. depend only on 
the deformation in its neighborhood. The corresponding 
short-range string has been the object of many theoret- 
ical and experimental studies. In the depinning limit, 
two different scenarios are possible; numerical simula- 
tions and analytical calculations [pUj] have established 
that a string with an elastic restoring force breaks at 
the depinning threshold, while percolation experiments 
and numerical studies on directed polymers in random 
media Mpl agree that in those systems with stronger 
than harmonic restoring forces the roughness exponent 
is C = 0.63. 

It has also been shown |||,|l^ that for a contact line 
of a liquid meniscus or for crack propagation in a solid, 
the elastic force is long-range, rather than local. Non- 
local elasticity can be expected to modify the dynamic 
and static properties of these systems and to change the 
critical exponents. In this work, we compute one of these 
exponents, the roughness exponent C, oia, long-range elas- 
tic string at the depinning threshold fc. 

The theoretical approaches are up to now based on the 
assumption that the motion of the line at the threshold is 
quasi-static. This means that velocity-dependent terms 
in the equations of motion of the deformation field h(x, t) 
are taken to be irrelevant and can be derived from an en- 
ergy function, which incorporates potential energy due to 
the driving force / and the disorder potential ri{x, h), as 
well as an elastic energy. According to this hypothesis, 
the equation of motion of the deformation field at zero 



temperature is: 

d 

— h{x,t) = f + T]{x,h) -fc 



h{x,t)-h{xi,t) 
^^1 1 ^2 • (1) 
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The last term in this equation accounts for long-range 
restoring forces. Let us note that measurements of lo- 
cal velocities for a liquid ^He contact line 0] have cast 
doubts on the validity of the quasi-static hypothesis for 
the present experiments. 

The critical behavior of eq.(|l|), at the driving force 
/ equal to fc, was studied by means of renormaliza- 
tion group (RG) techniques. The one-loop calculations 
|ll| gave a roughness exponent C}-^'^ equal to 1/3 which 
at a time was believed to be exact ||ll|,|l^. Early sim- 
ulations based on extremal long-range models [|l3|jl4[ 
(C = 0.35 ± 0.02) and on cellular automata H (( = 
0.34 ± 0.02) found good agreement with this conjecture. 
However, experiments, both in crack propagation |lq ] 
and for a liquid '^He contact line on a rough substrate 
y have measured, near fc, a systematically larger expo- 
nent (C — 0.56 ± 0.03). Chauve et al. recently showed 
that higher-order terms contribute to the RG result. At 
two-loop order, they found an exponent C*^^) = 0.47. 
The large difference of the two-loop calculation with the 
lowest-order result let it seem conceivable that the exper- 
imental value could be explained by the model eq. M) . In 
fact, the authors of ref. M estimated the exponent to be 
(^ — 0.5 ±0.1, which did include the experimental results. 

In this article we present a new and very precise 
method to determine the roughness exponent at fc. We 
show that the exponent ( has the value 0.390 ± 0.002. ^ 
is thus bigger than what was suggested by earlier simu- 
lations {cf however p7|) but is at the same time incom- 
patible with the experimental situation. 

The direct numerical simulation of the dynamics, es- 
pecially of long-range systems, in the depinning region 
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is extremely tedious because the velocity of the manifold 
vanishes at the threshold, which is thus difficult to ap- 
proach 1 18 19]. Many authors therefore preferred to treat 



the problem within the framework of cellular automaton 
models. These approaches can be very useful, but are 
hardly identifiable with a continuum equation (concern- 
ing the long-range case cf |l^,|2^ ) . 

However, there is an additional analytical structure in 



this problem, as first noticed by Middlcton |22| for the 
continuum models. For short-range lattice models |p|,pl[, 
we used this additional information to compute the crit- 
ical string (the blocked string at /c) in an extremely effi- 
cient way, but without actually simulating the time evo- 
lution of the system. On the lattice, the long-range model 
does not seem to be open to such an approach, but the 
continuum model is, as we will show in this paper. 
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FIG. 1. Example of a random potential on the i-th site 
with periodic boundary conditions (M = 32). The circles are 
Gaussian random numbers with zero mean and unit variance, 
assigned to evenly spaced points hi — 1,2,... ,M. These 
numbers are then interpolated by a periodic cubic spline, in 
order to yield a continuous random potential V{i, hi). 

To proceed, we discretize the variable x (x ^ Xi,i = 
1,... ,L) in eq.(|l|) and write h\ instead of h{xi,t). 
The complete string at time t is thus given by /i* = 
{h\,hl,. . . , /i^}, where the hj are real periodic variables 
for which hi+M and hi are identified. Periodic boundary 
conditions are also applied in L: hL+i = hi. The equa- 
tion of motion (Q) is adapted to accomodate the periodic 
boundary conditions. We have computed the long-range 
force f by summing over periodic images, as others have 
done before us |l^,|lj], obtaining: 



L-l 

E 

41=1 



f"'M = E 2^ («i/^) - vrVs*"'(^n/i) 



ft.* 



L2 



(2) 

where '0 = dT{x)/dx and T{x) is the gamma function. 
Different calculations, without the sum of images, gave 
an identical result for C. A random potential V{i,hi) 
acts on each site of the string. Our choice of the ran- 
dom potential is shown in fig.ll It allows to obtain a 
differentiable potential with r]{i,hi) — ~-dV{i,hi)/dhi. 



The 'no-passing' rule [g2[ establishes the following: if 
two strings h and h do not cross at a given time t (say, 
hi < hi Vi), they will not cross at any later time. An- 
other important property of Middleton's theorem states 
that if, at an initial time tinitj the velocities u| are non- 
negative for all i, they will remain so for all later times 
t > tinit- It follows from this property that, once we have 
found a forward moving string ft*>"'* ^ we can be sure that 
snapshots of the string at later times will never cross. In 
fact, the strings ft* for t > tinit will form a non-crossing 
family with non-negative velocities, which satisfies the 
following conditions: 



i) hi > hi' Vi for t > t' > tmit 

a) hi -^ hi' for t > t' > Unit, 

Hi) v{hl) > Vi for t > tinit 

iv) f^hl = v{hl) from eq.(§) 



t^t' 



(3) 



The velocity v(ft*) is given, as usual, by w(ft*) = 
—dE[h*]/dhl where Elh*^] is the energy of the config- 
uration ft*. In this context, we have made the follow- 
ing observation: if our only aim is to obtain the critical 
string, rather than to simulate the true time-behavior, it 
is sufficient to generate continuous non-crossing families 
satisfying i), ii), and Hi), without imposing iv). In this 
case, t would not be the physical time, but simply an 
ordering index; only then are the three conditions i), ii), 
and Hi) independent. 
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FIG. 2. a: 'No-passing' theorem: /i* (filled circles) and 
f^biock (-Qpgjj circles). If ft'* approaches h'''""'' in a point i, 
the disorder forces become equal, but the elastic term (repre- 
sented by arrows) prevents ft* from exceeding /i'''°'=* for t > t* . 
h: Algorithm. The point i of the string (filled circles) is 
moved from hi to hi in one step. Between hi and hi, the 
string's velocity in i remains positve and and is zero at ft^. A 
blocked configuration is drawn (open circles): our string may 
approach this configuration, but cannot pass it. 



To assure that we can drop the condition iv), we have 
to guarantee that no member of such a family will ever 
cross a blocked string ft^'°'='= (by definition v{hf°'^^) = 0), 
if ft*i"« did not. Let us suppose the contrary: if ft* were 
the ffist member which touches the blocked line in one 



point i (as shown in the fig.(ga)), the random force for 
j^biock g^jjj j^t ^Q^y become equal, but the elastic term 
would give w(/i* ) < v{h^^°^^). This would be a violation 
of condition in). Conversely, it is easy to see that such 
a family can always be continued up to a blocked line, 
because it suffices to find a single point i with positive 
velocity to continue the construction. 

In previous works pq |, the critical line was computed 
by direct simulation of cq.(|^). Notice that the discretiza- 
tion of time can pose difhcult problems: during the in- 
terval time Ai the motion does not respect Middlcton's 
theorem and therefore is not guaranteed to halt in front 
of a blocked string. 

The construction of continuous non-crossing families 
of strings presents a much more powerful, completely rig- 
orous approach. As indicated in fig.(Q), we consider in 
practice families in which at a given instant only a single 
coordinate i moves. Coordinate i is then advanced (from 
hi to h'j) until its velocity vanishes {v{h^) — 0). This is 
of course not done by simulation, but in a single step, 
by computing the zeros of v{hi). To simplify our numer- 
ical work, we have represented the continuous random 
potential as a spline (a piecewise third-order polynomial 
in /i), as further explained in fig.(n[). This allows us to 
solve for the closest zero of the velocity function from a 
quadratic equation. The only parameter in this proce- 
dure is the minimal velocity below which the total string 
is assumed blocked, and which serves to terminate our 
iteration. We have varied this non-essential parameter 
by four orders of magnitude, and shown that the critical 
forces and exponents are extremely well stabilized. 

An initial forward-moving string ft,*'"" is very easy to 
obtain. We have remarked that our iterative algorithm, 
which is completely rigorous, converges extremely quickly 
to a blocked string if such an object exists. A bisection 
method then allows us to obtain the critical driving force 
/c, and the critical string. We stress again that, by con- 
struction, a blocked string is never passed. 

We have run this algorithm on a large number of sam- 
ples with L = M. It can be shown that the correct 
thermodynamic limit is to consider change of scale as 
{L,M) -^ {aL,aM) pl|. Sample sizes varied from 
L = 64 to L = 1024, where we are able to obtain a 
critical string (at / = /c) in a few minutes on a PC. We 
compute the mean-square elongation W'^{L) of a critical 
string ft^ as 



W^{L) := < (ft<=- < h'' >)2 >. 



(4) 
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FIG. 3. Mean square elongation W'^{L) as a function of 
system size L for the long-range elastic force given by eq.(|l|). 
The interpolating line corresponds to a roughness exponent 
C = 0.390 ±0.002. 

In conclusion, we have obtained the roughness expo- 
nent for long-range elastic strings by a numerical proce- 
dure which respects the analytic structure in the problem 
('no-passing' theorem) and allows to obtain very high 
precision. Our approach can certainly be extended to 
other problems, we would in particular be interested in 
simulation methods which allowed to study, e.g., the KPZ 
equation ||23| directly in the continuum. This is most 
important, as it has been shown that the effect of dis- 
cretizations can be very difficult to control |Q . 

The large difference of the result with the experimental 
value shows that the theoretical model of these processes 
will certainly have to be modified in an essential way. 
One possibility is that velocity-dependent terms have to 
be taken into account. This appears reasonable, as large 
local velocities have been observed in the experiments. 
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In eq.(^, < h"^ >— j;J2i^iy while the overbar stands 
for an average over the disorder. Our data obviously 
extrapolate very well and a mean-square analysis yields 
C = 0.390 ±0.002. 
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